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Abstract 

We consider the intensity-based approach for the modehng of defauh times of one or more 
companies. In this approach the defauh times are defined as the jump times of a Cox process, 
which is a Poisson process conditional on the realization of its intensity. We assume that the 
intensity follows the Cox-Ingersoll-Ross model. This model allows one to calculate survival 
probabilities and prices of defaultable bonds explicitly. In this paper we assume that the 
Brownian motion, that drives the intensity, is not observed. Using filtering theory for point 
process observations, we are able to derive dynamics for the intensity and its moment generat- 
ing function, given the observations of the Cox process. A transformation of the dynamics of 
the conditional moment generating function allows us to solve the filtering problem, between 
the jumps of the Cox process, as well as at the jumps. Assuming that the initial distribution 
of the intensity is of the Gamma type, we obtain an explicit solution to the filtering problem 
for all t > 0. We conclude the paper with the observation that the resulting conditional 
moment generating function at time t, after Nt jumps, corresponds to a mixture of Nt + 1 
Gamma distributions. 

keywords: Credit Risk, Affine Model, Filtering, Point Process. 
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1 Introduction 



The main goal in credit risk is the modehng of the default time of a company or default times 
of several companies. The default times are often modeled using the so-called intensity-based 
approach as opposed to the firm value approach. Here, the default time of a company is modeled 
as the first jump time of a Cox process, of which the intensity is driven by some stochastic process, 
e.g. Brownian motion, or, in case of more than one company, as consecutive jump times of this 
Cox process. This approach enables one to calculate survival probabilities, and to price financial 
derivatives depending on the default of one or more companies, such as defaultable bonds and 
credit default swaps. The former pays a certain amount at the maturity of the contract in case 
the underlying company does not default, otherwise it pays a smaller amount, known as the 
recovery rate. The credit default swap is a form of default insurance, which pays the loss incurred 
on a default of the underlying company. Overviews of the intensity-based modeling approach can 
be found in [Lan98], [Gie04] and [EliOSa]. In this approach, it is a common assumption that the 
driving process can be observed, i.e. the observed filtration is generated by the Cox process, which 
can be seen as the default counting process, and by the driving process. 

In this paper it is assumed that the driving process is not observed, and thus only a point process 
Nt is observed, which introduces a stochastic filtering problem for point processes. In particular 
the intensity is assumed to follow the Cox-IngersoU-Ross (CIR) model, where the driving Brownian 
motion is not observed. General theory for filtering with point processes can be found in Bremaud, 
[BreSl], for example. The case in which one assumes a Cox-Ingersoll-Ross model for the intensity 
and that the initial value of the intensity is drawn from a Gamma distribution, has also been 
considered by Frey et al. [FPR07]. These authors derive a recursive solution to the filtering 
problem at jump times of the point process Nt. 

In contrast, in the present paper we also pay attention to the explicit solution to the filtering 
problem between jump times. We obtain this part of the solution analytically by solving partial 
differential equations. Furthermore, we consider a different approach to obtain the recursive 
solution at jump times. By combining these solutions, we obtain a solution for all t > Q. It is 
further observed that the resulting conditional moment generating function at time t corresponds 
to a mixture of A't -I- 1 Gamma distributions according to some discrete distribution. 
The paper is organized as follows: In Section 2 the Cox-IngersoU-Ross model is discussed and 
some results for the case of full information are discussed. Next, in Section 3, the filtering problem 
is introduced and some background is given for filtering of point process observations. First, 
the filtering formulas from [BreSl] are given, and the equations for the conditional intensity and 
conditional moment generating function are derived. Then, in the second part of Section 3, we 
introduce filtering by the method of the probability of reference, and the filtering equations are 
transformed using the ideas introduced in [BB80]. Section 4 deals with the filtering problem 
between the jump times of the point process, given the initial distribution of the intensity at 
jump times. In Section 5, the filtering problem is solved at jump times, and an explicit, recursive 
solution is obtained, which combines the solutions between and at jumps. Further the resulting 
conditional moment generating function is analyzed and it is observed that this function agrees 
with the moment generating function of a mixture of Gamma distributions. The section concludes 
with an illustration of the mixing probabilities. 

2 Model and Background 

One of the main goals in credit risk is the modeling of the default time of a company or the default 
times of several companies. Over the years two approaches have become popular, the structural 
approach and the intensity-based approach. In the structural approach the company value is 
modeled, for example as a (jump-) diffusion, and the company defaults when its value drops below 
a certain level. This approach is discussed in more detail in e.g. [Gie04], [BR02] and [Eli05b]. In 
the intensity-based approach the default time is modeled as the first jump of a point process, e.g. 
a Poisson process or, more general, a Cox process, which is an inhomogeneous Poisson process 
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conditional on the realization of its intensity. In case one considers more than one company, one 
can model the default times as consecutive jump times of the Cox Process. In [Lan98], [Gie04] 
and [EliOSa] this modeling approach is discussed in more detail, and [Sch02] provides a detailed 
application. In this paper we focus on the intensity-based approach, where the intensity, Xt, of 
the Cox process has an affine structure, similar to interest term structure models [DK96]. This 
means that the intensity process At follows a stochastic differential equation (SDE) of the form: 



dXt = (a + bXt)dt +y/c + dXt dWt, (2.1) 

for a Brownian Motion Wt, with > 0. In particular, the focus is on the Cox-IngersoU-Ross square 
root (CIR) model, [CIR85], for the intensity, where the intensity At satisfies 



dAt = -a(At -Aio)dt + /3VAtdWt. (2.2) 

In [LL96, Section 6.2.2.] one finds parameter restrictions for this model which guarantee positivity 
of At- Naturally one should start with a positive initial value Aq, and if afiQ > /3^/2, then Af 
remains positive with probability one. Note that using the transformation Xt = Xt+c/d and by a 
reparametrization, Xt satisfies the general SDE (2.1), and At satisfies (2.2). This implies that the 
general form (2.1) and the CIR intensity (2.2) arc in fact equivalent. Therefore the CIR intensity 
will be considered in the remainder of this paper. 

A big advantage of the afhne setup is that many relevant quantities in credit risk can be calculated 
explicitly. Using the formulas from [LL96, Section 6.2.2.] one can, for example, easily calculate 
the survival probability P(t > t\!Fs), with t > s and !Ft = V J^j^, where the former filtration is 
generated by the point process Nt and the latter by some process Yf driving the intensity process. 

Example 2.1. Consider, on the filtered probability space {fl, {Tt)t>o, a random time r > 
as the first jump time of a Cox process Nt, which intensity follows the CIR model (2.2). Further 
assume that .Ft = V , where J-^ is the filtration generated by the Brownian motion that 
drives the intensity process. Then one can calculate the survival probability ior t > s as 



(r > t\J^s) = l{r>s}]l 



(2.3) 



which follows from formulas in [BR02, Chapter 6]. Since At is a Markov process, one can condition 
on As instead of J-^. An application of Proposition 6.2.4. from [LL96] to (2.3) yields 

P(r > t\Ts) = l{r>s} cxp i-afio<P{t - s) - X,ij{t - s)) , (2.4) 

where 

0(0 = -^ log' 



^2 ^ ^ a + e*"' {'^ + a) 
2 (e"^* - 1) 



7 — a + e*'>'(7 + a) 



7= 

Other relevant quantities, such as the price of a defaultable bond, can also be calculated analyt- 
ically, under some restrictions on the interest rate, e.g. by posing that the interest rate evolves 
deterministically. In [FPR07] some of these quantities are considered in more detail. 
It is a common assumption, which is also followed above, that the filtration J-t is built up using 
two filtrations, and , where the first filtration represents the information about the process 
driving the intensity and the second filtration contains information about past defaults. In this 
paper it is assumed that the factor Y is not observed which results in a filtering problem of a point 
process. 

In the following sections the problem is introduced formally and solved for the case where the 
intensity follows the CIR model. 
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3 The Filtering Problem 



In filtering theory one deals with the problem of partial observations. Suppose that a process Zt 
on the probability space P) is adapted to the filtration Tt- Furthermore let the process 

Yt be observed, where Yt is measurable with respect to a smaller filtration Tl[ C Tt- One is 
then interested in conditional expectations of the form = E \^Zt\J-^], and one tries to find 
the dynamics of the process Zt, for instance by showing that it is the solution of a stochastic 
differential equation. 

In this section the filtering problem is considered in the case a point process is observed. First 
some general theory about filtering with point process observations is discussed, and Example 2.1 
is continued within the filtering setup. The calculation of the survival probability depends on the 
conditional moment generating function, for which an SDE is derived. In the second part of this 
section this equation is transformed in such a way that the filtering problem allows an explicit 
solution. 

3.1 Filtering Using Point Process Observations 

In the case of point process observations the observed process Yj is equal to the point process Nt, 
with ^(-intensity Xt. The process Zt is assumed to follow the SDE 

dZt ^ atdt + dMt, (3.1) 

for an J-'t -progressive measurable at, with |as|d.s < oo, and an J-"t-local martingale Mt- The 
filtering problem is often cast as the calculation of the conditional expectation E[Z(|jFj^] =: Zt. 
Using the filtering formulas from [BreSl, Chapter IV], a representation of the solution to this 
filtering problem can be found. In case the (local) martingale Mt and the observed point process 
have no jumps in common, one has: 

dZt = atdt + (^^Y^ " ^*-) i'^^' ' ^"^*) ' ^^-^^ 
with at ■■= ElatlT^], and Xt- := liuis-it Xg. 

Example 3.1 (Example 2.1 continued). When one wants to calculate the survival probability 
given T^ , one has Zt = l^T>t}- Combining this with the survival probability in the case of full 
information, one can calculate the survival probability P(r > t\T^). 

P (r > t\Tf) = E [P (r > t| J-f V T^') I ] 

= l{r>s} exp(-Q;^o0(i - s))E [exp{-ip{t - s)As)| J^f ] , 

which can be calculated if an expression for the conditional moment generation function f{s, t) :~ 
E [e"^*| J^f"] is available. 

The above example illustrates that one can calculate the survival probability if the conditional 
moment generating function /(s, t) is known. As a first step in the determination of this function, 
the SDEs of At := E [At| and f{s,t) are determined. First Ito's formula is used to obtain the 
SDE for e"^' , where A* satisfies (2.2) 



-as + Is^pA ^e"^' + sa^ioe'^* 
2 J as 



dt + |3^/Xte'^'dWt. 



The filtered versions are obtained by applying formula (3.2). One obtains for A 



dAt = -a(Xt - Mo)dt + f ^ _ At_ j (^dNt - Xtdt ) , (3.3) 
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and for /(s,t) one finds 
1 



df{s,t) 



d 



^ ) -g^fis,t) + saiiofis,t) 




Atdt 
(3.4) 



In general, filtering equations are very difficult, if possible at all, to solve explicitly, since the first 
equation involves terms with and the second equation involves combinations of At and f{s,t). 
In order to solve these equations one should also have equations for A^ , but this involves A^ and 
so on, assuming that they exist. So instead of trying to solve these equations directly, a different 
approach is considered in order to find an expression for /(s, t). 

3.2 Filtering by the Method of Probabihty of Reference 

In order to solve the problem introduced above, the filtering by the method of probability of reference 
is considered, see [BreSl, chapter VI] or [BB80, Section 2]. In this approach a second probability 
measure Pq and intensity process A^ are introduced, such that Nt — /J X^ds is a martingale with 
respect to J-'t under Pq. Corresponding to this change of measure one has the likelihood ratio, or 
density process A, given by 



At :=E 



" dp 












[dPo 


Jo 



A,_ 



A? 



— (dNs - X"ds) . 



(3.5) 



This likelihood ratio turns out to be a useful tool to solve the filtering problem for f{s,t). It 
is known, see e.g. [BreSl] for the case Aj = 1, that the filtered version of this likelihood ratio. 
At := E [At| ^[^] follows an equation similar to (3.5). One has 



At = 1 



A. 



As--Ag_ 
A9 



diV, - A' 



To solve the filtering problem for f{s,t) an auxiliary function g{s,t) is introduced. It is defined 

by 

g{s, t) f{s, t)At exp ( - / A°dM 



(3.6) 

The exponent is used in order to obtain a simpler SDE of g{s, t). After a solution to this equation 
has been found, one can obtain f{s,t) by 



f{s,t) 



5(0,0' 



(3.7) 



It is directly clear that the first and third component of g{s,t) are positive, and from (3.9) follows 
that also the second component is positive, and thus the division in (3.7) is well defined. The 
solution to the filtering problem is obtained as soon as an expression for g(s, t) is found. In 
Proposition 3.2 an SDE is derived for g{s,t) for the intensity following the CIR model. 

Proposition 3.2. Let g{s,t) be given by (3.6), then one has, fort > 



dg{s,t) 



sfiQag{s, t) + ( ^s^P^ 



dt 



'1 _a 
ds 



9is,t-) - g{s,t-) 



dNt. 
(3.8) 



Proof. As a first step in proving (3.8), one can rewrite the function g{s,t). An alternative expres- 
sion for At is given by 



At 




cxp 



T„<i 



A„ 



A" ) du 



(3.9) 
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which can be checked by a direct calculation. From this it is easy to see that 



g{s,t) ' 7(s,i)AtCxp (^-^ Zdu^ = /(.s,t) J] 



exp - / Xudu =: f{s,t)Lt 







dLt = ( dNt - A?dt) - Lt-dNt. 



For Lt one finds the SDE 

Lt-Xt- 

The SDE in (3.8) foUows from the product rule 
d.9(.s, t) = /(.s, t-)dLt + Lt-df{s, t) + A/(s, t)ALt 



: f{s, t-) (^^^^^ (dA^t - ^tdt) - Lt-dNt^ + Lt- (^-as + is^/?^^ ^/(s, t) + sa/^ofis, t) 

+ ( - /(. t-)] (dNt - xtdt) ] + ( - /(. f - 2.-1 d^. 



dt 



Collecting the terms before dt and dA''^, one obtains the equation 

dg{s, t)= ^tfi-s, t)Lt + {-as + is^/?^^ ^/(s, t)Lt + saAio/(s, t)Lt - ^J{s, t)Lt + f{s, t)Lt\}j dt 



\t- Xt_ 

d _ 

+ f{s,t-)Lt- ]dNt. 



f(s,t-)Lt-Xt- §-J{s,t-)Lt 



A?_ Xt- 

The result follows by simplifying the last equation. □ 

The right hand side of (3.8) depends only on g{s,t) and its partial derivative with respect to s. 
In the next section this equation is solved between jumps, and in section 5 the equation is solved 
at jump times of the process Nt- 

4 Filtering Between Jumps 

In the previous sections the filtering problem for point processes has been defined in general 
terms, and the problem has further been considered for an intensity following the Cox-IngersoU- 
Ross model. To solve the filtering problem, one has to solve equation (3.8). This equation can be 
split up into a partial differential equation between jumps of the process Nt and an equation at 
jumps. In this section the equation between jumps is solved for a general initial condition at time 
T > 0. Later on T will be considered as a jump time of Nt- Note that an initial condition for 
g{s,t) is given as 

g{s, T) - /(s, T)Kt cxpl- [ A^du ) . 



For T = it follows that 

5(s,0) = /(s,0) = E [e'^°\T^] E [e^^«] , 

which is the moment generating function of the intensity at time t = 0, since J-^ — {0, il}. 
Before the solution to (3.8) is found, the specific case is considered, in which all the parameters 
in the CIR model are set to zero. Albeit a simple example, the analysis of it sheds some light on 
the approach that will be followed for the general case. 
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Example 4.1. Consider the CIR model in which all the parameters are set to zero. This results 
in a constant intensity, and thus dXt = 0. The filter equations (3.3) and (3.4) reduce to 

dXt= (^j-^Xt-^ (dNt ~ Xtdt 

dft = (^^y^ - ft-^ {dNt - Xtdt 
The partial differential equation for g{s, t) between jumps reduces to: 

^^g{s,t) = -^^g{s,t). 

With an initial condition g{s, T) ~ w{s), one easily finds that the solution to this equation is 

g{s,t)^w{s-t + T). 

In the next section this example is considered once more, where the filter at jump times is consid- 
ered. We proceed with the case of an intensity following the CIR model. 

Proposition 4.2. Let At follow the Cox-Ingersoll-Ross model (2.2), and let g{s,t) be given by 
(3.6), with an initial condition at time T , g{s,T) = w{s). Then, for T <t < T^, with T„ the first 
jump time of Nt after T, g{s,t) solves the partial differential equation 

d 1 d 

— 5(s, t) ^ sfioag{s, t) + —{ps - a + T){ps - a - T)—g{s, t), (4.1) 
ot 2p OS 



where p := P"^ and r :~ -^/a^ + 2/3^. The unique solution to this equation is given by 

n(s t) = e«("-^)(*-^) ( — 

^ps(e--(*-^)-l) + (T-a)e— (*-^)+T + ay 

/s ((a + r)e--(*-^) + T - a) + 2e--(*-^) - 2\ 
ps(e— - 1) + (r - a)e— (*-^) +T + a ) ' ^^'^^ 

where 61 

p 

Proof. The partial differential equation (4.2) for g{s, t) follows directly from Proposition 3.2, since 
the jump part of this equation can be discarded. 

To obtain a solution to this equation a candidate solution is derived by making a number of 
transformations of the independent variables, until a simple PDE is found, which can be solved 
explicitly using known techniques. This candidate solution can then be checked to be the solution 
by calculating its partial derivatives, and inserting these into (4.1). 
The first transformation is given by 

{s,t) (Pl^^±l,t] {u,t). (4.3) 
\ps — a — T J 

Instead of g{s,t) one writes /i(u, s), in terms of the new variable u. Using this transformation 
and the PDE for g{s,t), one can derive a PDE for fi{u,t), by expressing s in terms of u, and 
expressing the partial derivatives of g{s, t) as partial derivatives of fi{u, t). The resulting PDE for 
fi{u,t) is 

9 . , , fa t(u + 1) \ ^ . . d . . , 

-f,{u,t) ^ MOO + ) f.M - ru-f,M. 

The second transformation that is used is given by 
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where, for the thne being, u is tacitly understood to be positive. Instead of the function /i(s,t), 
one considers the function f2{v,t) := /i(u,<), in terms of the new variable v. This transformation 
results in a partial differential equation for f2{v, t), 

The final transformation is given by 
which results in the PDE for /3(w,t): 

This equation can be solved using the method of characteristics, which is explained in chapter 1 
and 8 of [Che71], for example. Using this technique the partial differential equation is transformed 
in an ordinary differential equation by introducing new variables £,{v,t) and The former 

is used to replace both v and t, and the latter is used to parameterize the initial curve. To be 
able to solve the PDE an initial condition is required for /^{vjt). By applying all the previous 
transformations to the initial condition g(s, T) = w{s), with t >T, one obtains the initial condition 

Next one has to solve the differential equations 

with the initial conditions t(0, Q) = T and w(0, Q) = Q. The unique solution to these equations is 
trivially given by 

Inverting these expressions, yields 

av,t)^t-T, Civ, t)=v-t + T. 

Using these transformations, the partial differential equation (4.4) can be transformed into the 
ordinary differential equation (ODE) 

:/3(e,C)=Moa - + -7-Z77Tn 7T = : + -J-^TTTn 7T-^(" + ^) + - 



(4-5) 

where 6 = H^Jl^ This ordinary differential equation can be solved for the given initial condition 
f^(v,T) = G{v). To derive the solution one can start with a candidate solution 



/3(e, C) - Ci log (e"(«+';) - l) + CaC + C3. 



For ^ = 0, one has /3(0, C) = C*! log (e"^^ — l) + C3, and /s has partial derivative with respect to 

^/3(e,C)=rC, + ^,(^^^ + C,. 
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Using the initial condition /3(0, C) — G{C), together with the ODE (4.5), one can find the values 
of Ci, C2 and C3: 



Ci = 20, 

C2 = 0{a - r), 

C3=G(C)-201og (e-«-l). 

This leads to the unique solution 

/3(e,C) =e(a-T)e + 201og(e^(«+^)-l) +G(C)-201og(e^«-l). (4.6) 

The proof of the uniqueness of this solution is postponed to the end of this proof. 
Replacing ^ by i — T and C by u — < + T, results in 

Mv, t) = 0(a-r)(t-r)+2^?log ^^,1^^^, ^ +log w _,JjJ \, • (4.7) 



Next, one obtains a candidate solution for 17(5, t), by reversing all the transformations. This gives 



u-1 /ue-^(*-^)(r + a) + r 



By performing the last substitution, (4.3), an expression for g{s,t) is obtained. One has 



/ ps-a+T _ -, \2e / ps-a+T rjt-T}/ , , . 



ps — a — T 



gS(a-T)(t-T) 



\ps — a — T J 

2r 



X w 



Vps(e— ^(*-^) - 1) + (r - Q;)e— ^(*-^) + T + a 

's{{a + T)e-^(*-^) + T - a) + 2e-(*-^) - 2\ 
ps(e-^(t~T) _ 1) + _ Q,)e-^(*-^) +T + a I ' 



where it was used that (a + t){t — a) = 2p. By inserting this candidate into equation (4.1), one 
can check that it indeed is the solution. 

The last thing to proof is the uniqueness of the solution to equation (4.1). As all the transfor- 
mations are clearly one-to-one, the uniqueness of this solution should follow from the uniqueness 
of the solution to equation (4.5). It is easy to see that the solution to this equation is imique, as 
the difference of two possible solutions, with the same initial condition, has zero derivative, which 
implies that the two solutions are in fact equal. □ 

The result of Proposition 4.2 tells us that one can calculate g{s, t), for T < t < Tn, where r„ is the 
first jump time of Nt after T. In order to completely solve the filtering problem, one further has 
to solve the equation (3.8) at jump times. This is the topic of the next section, where a recursive 
solution will be obtained for the case in which Aq has a Gamma distribution. 



5 Filtering at Jump Times and a General Solution 

In the previous section the filtering problem has been solved between jumps, for an arbitrary 
initial condition w{s) for g{s,t), at time T > 0. In this section the filtering problem is solved at 
jump times, first for Example 4.1, and after that for the case where the intensity follows the CIR 
model. 
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Example 5.1 (Example 4.1 (continued)). At jumps one obtains from equation (3.8) 



Agi-s, t) = y^^:^^^ - gis, t-) j AN^. 

From this identity it easily follows that at a jump time T > 0: 

gis,T)={\'^_y'-^^gis,T-). (5.1) 
Combining the results between jumps and at jumps, one can obtain the solution to equation 

dg{s,t) ^ ^^g{s,t)dt+ ( - g(s,t-)] dNt. 



t- 



At each jump time T„, one has to take the derivative of the function g(s,t), and divide by Ai^ 
the resulting function can then be used as initial condition for the interval [T„,T„+i). Using an 
initial condition g{s, 0) = w{s), one obtains the solution 



5(s,t) = ^«(^')(s-i)n A 



n=l 



^0 



where w*-"^(s) denotes the rt-th derivative of w{s). The conditional moment generating hmction 
is found from (3.7), and is given by 

If one assumes that Aq r(a, /?), one has 

/(.,o) = 7(.,o)= (^)", /(^'^)= (^^)""- (^-2) 

From this follows that at time t > 0, Xt given is distributed according to r(a + Nt,(3 + t). 
Further At can easily be derived by a differentiation with respect to s: 



OS 



_ a + Nt 



The solution in this example was easy to find, which could be expected, since At is constant over 
time in this case. The general Cox-Ingersoll-Ross model for the intensity is more complicated, but 
in the remainder of this section, also this problem is solved. At jumps one has the same equation 
as in Example 5.1, which is already solved in (5.1). In Theorem 5.2 the solution for g{s,t) for the 
CIR model is given. Before this theorem is stated some notation is introduced. 
Let x,y & M., t > and put 

A{x, t, y):=x{{T~ a)e-^* + r + a) + 2y (l - e"^*) (5.3) 
B{s, t) ps (e^^* - 1) + (r - a)e"^* + r + a (5.4) 
C(x, t, y) y {{a + T)e-^* + t - a) + px {l - e"^*) . (5.5) 

This notation allows us to write the general solution between jumps, (4.2), as 
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Next let Ti, T2, . . . denote the jump times, and let Tq — 0. Then introduce the following notation: 



^(t, To) 



A{(l>,t,l) iorQ<t<Ti, (5.7) 

A(^(T„,T„_i),i-T„,^(T„,T„_i)) forT„<i<T„+i, (5.8) 

C(0,t,l) forO<i<Ti, (5.9) 

c(^(r„,r„_i),t-r„,^(r„,r„_i)) for t„ < i < t„+i. (5.10) 



With this notation, the main result of this paper can be stated. A recursive solution to the filtering 
problem is obtained, for the case where Ag has a Gamma distribution. 



Theorem 5.2. Let Xq - r{29, (p), for (p > and 9 = > 0. Then one has 

fo{s) = g{s,0) -- 



p 

28 



— S 



which is the moment generating function of the T(26, (j)) distribution. With the notation introduced 
in (5.3)-(5.5) and (5.7)-(5. 10) one further has, for Tn < t < Tn-^-i, 



^(-'^)-^(^)^"(-'n ^(^,T„)-.^(t,r„) j ' ^'-''^ 

where po{s,t) = 1, and for n > 1, pn{s,t) is a polynomial of degree n in s, that satisfies the 
recursion. 



c(-|,t-T„,s) 



Pn{s,t) = B"{s,t-Tn) 

(5.12) 

where di denotes the derivative with respect to the first argument of Pn, and 

'(2r)2e^ 

m>l, T„<t \ ^T,„- 

In the proof of this theorem the following lemma is used. 



X(t) = e«(— )*(2r0)^^ n (5-13) 



Lemma 5.3. With the notation from (5.3)-(5.5) and (5.7)-(5.10), the following relations hold for 
n > 1 and x,y Cz M; 

(i) ^/(T„,T„) = 2r.e/(r„,r„_i) 
(ti) '^(T„,T„) =2t<^(T„,T„_i) 

(Hi) xB{s,t) -yC {^-j^,t,s^ ^ A{x,t,y) - sC{x,t,y). 
Proof. (i) From equations (5.8) and (5.3) follows that 

^(r„,r„) - A(^(r„,r„_i),o,'^(T„,T„_i)) 

- ^/(r„, r„_i) ((r ~ a)e" + r + a) + '^^(T„, T„_i) (l - e") 
= 2r^(T„,T„_i). 

(ii) This follows along the same lines as in (i), using equations (5.10) and (5.5). 
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(iii) Using equations (5.3), (5.4) and (5.5) one finds: 

xB{s, t) - yC ( --,t,s] ^xips (e^^* - l) + (r - a)e"^* + t + a) 



P 



-y{s {{a + r)e-^* + r - a) + 2 (l - e"^*)) 
--x{{t- a)e~^* + T + a) +2y{l- e"^*) 

- s{y {{a + r)e"^* + T-a) +xp{l- e"^*)) 
^ ^(a;,t,2/) - sC{x,t,y). 



□ 



Now, Theorem 5.2 can be proved. 



Proof of Theorem 5.2. For each n it has to be shown that (5.11) holds at r„, and between r„ and 
T„_|_i. First this is shown for n ~ 0. Then the induction step is proved for n > 1. 
n = 0: For t = Tq = one has by assumption: 



From (5.11) one finds: 

5(s,0) = A'(0)po(s,0) 
= e"(2T0)2« 



2 \ 28 



^(0,0) - s<^(0,0) 



2 \ 28 



/ 2r(? 



^(0,O,l)~sC(0,O,l) 

28 / J X 28 



\2T(i)-2TS 



Next the interval up to the first jump time, < <: < Ti, is considered. From (5.6) and the 
expression for w{s) — g{s,0), one finds: 



28 




Kit)pois,t) 



^ X 28 



which is the same expression as in (5.11) for n ~ 0. The final step in the derivation above follows 
from Lemma 5.3 (iii), with x — (j) and y — 1, together with the definition of K{t) in (5.13). 

n > 1: Now it remains to prove the induction step. Therefore one can assume that equation (5.11) 
holds for n — 1. It then remains to show that the equation holds for n, at r„ and between r„ and 
T„+i. First the jump is considered. Thus one has to calculate the derivative of 17(5, t) with respect 
to s, and take the left limit in t = Tn, further the derivative is divided hy _. By (5.1) one has 

g{s,Tn) = [^T„-) -g^9is,Tn-) 

10/ / 1 \ 28+n-lN 



£/(r„,r„_i) -s^(T„,T„_i) 
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Calculating the derivative with respect to s, leads to 



p„_i(s,T„)(20 + n-l)'r(T„,T„_i) 



— p„_i(s, r„) (^(T„, T„_i) - s<^(r„, r„_i)) 



^(T„,T„_i)-s'^^(T„,T„_i) 



2e+n 



(5.14) 



From Lemma 5.3 (i) and (ii) follows that for the denominator in (5.14) one has 

^(T„, T„_i) - s'^(T„, r„_i) = (2r)-i (^(r„, T„) - s'^(T„, r„)) . 
Hence (5.14) can be written as 



5(s,T„) = /Y(r„-)(2T)2''(2r) 



+ — p„_i(s,T„) (£/(T„,T„_i) ~ .s<r(T„,T„_i)) 



p„_i(s,T„)(20 + n - l)'r(r„,T„_i) 

1 



£/(r„,T„)-s'^(r„,T„) 



26+71 



(5.15) 



From (5.13) it is easy to see that A'(r„) = A'(T„-) (^t„-) (2t)^'', and frn'ther one has 2t = 
i3(s, 0) = i3(s, Tn — Tn). From this follows that (5.15) can be written as 



5(s,T„)=i^(T„)B"(s,r„-r„) 



Pn-lis, T„){2e + n- 1)^^(T„, r„_i) 



+ — p„_i(s,T„) (£/(T„,T„_i) - .s^(T„,T„_i)) 
OS 



^(r„,r„)-s'^(r„,T„) 



This can be simplified further using the definition of p„(s,t) as given in (5.12). together with the 
identity C f — ^,0,sj = rs. This results in 



g{s,Tn) = KiTn)p„{s,Tn) 



^(r„,r„) -s<r(T„,r„) 



28+71 



which is the required result at t = Tn- Finally one has to check that (5.11) holds for Tn <t < Tn+i- 
For this one can use the general solution (5.6) with initial condition w{s) = g(s, Tn). One finds 



g{s,t) = e 



e(a-T)(t-T„) 



2r 



B{s, t - Tn] 



2B 

ge(a-r)T„^2r</))2e [] 



(2r 



,20 



S(s,<-T„ 
Simplifying this expression yields: 



1 



m>l, T,„<T„ \ '~T, 



1 



20+n 



^S^(Tn^Tn) ^B{s t—Tn) ^ i'^fi j ^7i) ^ 



g{s,t) = e<'i--^^\2rc^r R (|r^jj (2r)^''S"(., t - T„)p„ 



I t — Tn, S 

B{s,t-Tn) 

26+71 



2tB(s, i - T„)i/(T„, T„_i) - 2rC ( -|, t - r„, s ) '^(r„, T„_i) 
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An application of Lemma 5.3, with x = ^(r„,T„_i) and y = "^(Tn, T„_i), and the definitions of 
t,Tn) and '^{t,Tn) in (5.8) and (5.10), together with the definition of K{t) results in 



Next, with the definition of p„(s, r„) from (5.12), evaluated in t = T„, together with C{x, 0, y) 
2Ty and B{s, 0) = 2t one rewrites this to 



g{s,t) = K{t) 



1 



Pn-1 




B{s, t - T„) 

B{s,t-Tn) 
1 

^(i,T„)-s'^(<,T„) 




'^(T„, Tn-l) 



In the final step the definition of p„(s,t) is used, this time evaluated in i, which concludes the 
proof of (5.11). From the definition of B{s,t) and C{x,t,y), with y = s, which are both linear in 
s, follows that Pn{s, t) is a polynomial of degree ri in s. □ 

This theorem provides a recursive solution to equation (3.8), in case Aq is distributed according 
to a r(20, 0) distribution. From (3.7) it already known that the conditional moment generating 
function can easily be obtained from an expression for g(s,t). Now this has been found, the 
conditional moment generating hmction /(s,t) can be obtained easily. 

Corollary 5.4. Under the assumptions of Theorem 5.2 the conditional moment generating func- 
tion f{s, t), for Tn < t < Tn+i, can be expressed as: 

f{s,t)^q^{s,t) ( y^y (5.16) 



\Q{t,Tn)-S 



where 



Here qn{s, t) is a polynomial of degree n in s. 

Proof. The result follows directly from equation (3.7), Theorem 5.2 and the definitions of g„ and 

g. □ 

With the derivation of the conditional moment generating function the filtering problem has 
been solved, and one is able to calculate conditional default probabilities using the results in 
Example 3.1. To conclude this section it is observed that the conditional moment generating 
function in (5.16) corresponds to a mixture of Gamma distributions. 

Remark 5.5. Corollary 5.4 provides an expression for /(s, t) that involves the polynomial (?„(■, t). 
Deriving an explicit expression for qn{s,t) = pn{s,t)/p„{0,t) for any n > is quite complicated, 
but we can write 



.Xs,t) = j2Rnt)s\ 
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where the coefficients, Rfit), of the polynomial follow directly from the coefficients of the polyno- 
mial in s, Pn{s,t), which in turn can be obtained using the recursion (5.12). 

Next, one can consider n + 1 independent random variables Ti, where Ti ~ r{2d + n — i,Q{t, T„)), 
for i ~ 0,1, ... ,n. Further, consider the discrete random variable A/", independent of the Ti, 
which assumes the values 0,1, ... ,n, with probabilities tt" (t), and define the random variable 

n 
1=0 

The moment generating function of X" can easily be found, as Ti and M" are independent, hence 

,Q{t,Tn) - 



E 



e-"| =EEK^'l{M=d ^Y.<m-''^] -E<W (tJ^) ■ (5-17) 



4=0 i=0 4=0 



The goal is to show that, by choosing the probabilities correctly, the moment generating function 
of X" equals the conditional moment generation function f{s,t). Therefore (5.17) is first rewritten 
as 



E 



/ Q{t,T^) " / Q(t,r„)-. 

{Q{t,T,,)-s) £^ Q{t,T,,) 



To have that both moment generating hmctions f{s,t) and (5.17) are equal, it is required that 

i=0 i=0 ^ ^\ y nj 

The right hand side of this equation can be written as 

i=0 j=0 ^ ■' ^ 

This equation can be turned into a polynomial in s, by interchanging the summations, which leads 
to 

n n y . \ n I n / . 

EE ( j) ^7m{t.T,.)-'^'{-iy ^Y.^' {-iyQ{t,T.n)-^Y.( ] ) ^rw 

j=0 i=j ^ ^ j=0 y i=j 

The moment generating functions are equal when 

R-{t) = {-iyQ{t,T,-,)-^j2(] )7r,r(i), 

for i ~ 0,1,..., rt. This can be solved iteratively, starting from j — n, which results in the 
probabilities 

n^lit) = {-iyR-{t)Q{t,T,,y - {])■ (5.18) 

\ ■' / 

It is not immediately clear from (5.18) that the 7r"(t) are all non-negative and sum to one. It 
turns out however that this is indeed the case for Tn < t < Tn+i, which means that the 7r"(t) 
can be interpreted as probabilities. It is however far from trivial to provide a general proof for all 
n > 0. We confine ourselves to illustrate this fact by some examples. In figure 1, two graphs are 
given in which the probabilities arc plotted. 
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1.0 r 1-0 r 




34567 45678 

(a) tTq (short dashed line), (long dashed line) (b) tTq (dot dashed line), Trf (short dashed line), ttj 
and 7r| (solid line) (long dashed line) and 7r| (solid line) 



Figure 1: Graphs of the mixing probabilities after two jumps of the process Nt, (a), and after three jumps, (b). 
The values of the previous jump times, Ti and T2 in case (a), and Ti, T2 and T3 in case (b), are taken as Ti = i, 
such that one is able to calculate the ■n"{t). The model parameters are chosen to be « = 0.5, f) = 0.5, ^0 = 0.4 
and 4) = 4.0. 
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